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)£ • Abstract 
Q_i| The scalar potential of the MSSM may have local and global minima characterized by 



non-zero expectation values of charged and colored bosons. Even if the true vacuum is 
not color and charge conserving, the early Universe is likely to occupy the minimum of 
the potential in which only the neutral Higgs fields have non-zero vev's. The stability of 
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this false vacuum with respect to quantum tunneling imposes important constraints on the 
values of the MSSM parameters. We analyze these constraints using some novel methods 
for calculating the false vacuum decay rate. Some regions of the MSSM parameter space 
are ruled out because the lifetime of the corresponding physically acceptable false vacuum 
is small in comparison to the present age of the Universe. However, there is a significant 
fraction of the parameter space that is consistent with the hypothesis that the Universe rests 
in the false vacuum that is stable on a cosmological time scale. 
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1 Introduction 



In the Standard Model color and electric charge are automatically conserved because the 
only fundamental scalar field is the Higgs boson, a colorless electroweak doublet. The Higgs 
potential has a continuum of degenerate minima, but these are all physically equivalent, 
and without loss of generality one can always define the unbroken U(l) generator to be 
the electric charge. This is not the case in the Minimal Supersymmetric Standard Model 
(MSSM), which employs a pair of Higgs doublets as well as a number of other scalar fields, 
the supersymmetric partners of quarks and leptons. Although the relative alignment of the 
two Higgs doublets' vacuum expectation values in group space is physical, the minimum of 
the Higgs potential (at least at tree level) preserves electric charge[| as long as the squark 
and slepton fields have vanishing classical values (see, e. g., Ref. JT] and references therein). 
However, the full scalar potential of the MSSM may have additional charge and/or color 
breaking (CCB) minima due to the vacuum expectation values of charged and/or colored 
scalars. 

The existence of the CCB minima in the MSSM in addition to the acceptable standard- 
model-like (SML) minimum may have important physical consequences. One might expect 
that the regions of parameter space for which there is a global CCB minimum could be 
automatically excluded, thereby further restricting theoretical predictions for the MSSM 
spectrum. However, one must be careful in drawing such conclusions. Just as the cup being 
the lowest point on the golf course by no means guarantees that the ball will end up there 
after being struck, the Universe at present may not be in its lowest possible energy state. 
Instead, it may rest in a false vacuum whose lifetime is large on a cosmological time scale. The 
fundamental reason that makes this possible is that quantum tunneling, a non-perturbative 
effect responsible for the first-order phase transitions in field theory, naturally introduces a 
time scale that is exponentially larger than the typical scale that characterizes the effective 
potential. Consequently, the relaxation to the lowest energy state from some excited state 
may take a very long time. In particular, parameters for which the local SML "false vacuum" 
1 This need not be the case in a general model with two Higgs doublets. 
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has a lifetime large in comparison to the age of the Universe may be acceptable, provided of 
course that the SML minimum was populated first in the evolution of the Universe. 

The existence of local CCB minima which were populated temporarily during the early 
stages in the evolution of the Universe would also have dramatic implications for cosmology 
and astrophysics. In particular, since baryon and lepton numbers are spontaneously violated 
in the CCB vacua, their existence might have important consequences for baryogenesis. 

Previous attempts 0-0 to elucidate the structure of the CCB minima in the MSSM 
met with serious difficulties. Some analyses 0, §], |5|, §] attempted to find analytic constraints 
on CCB minima. However, such conditions are generally neither necessary nor sufficient 
0] except for overly simplified toy models which resemble the MSSM in some features, but 
cannot be used to draw firm conclusions about the MSSM. For this reason, recent studies 
H 0, |] have employed extensive numerical analyses. Secondly, the determination of whether 
a global minimum is "dangerous", or not, must rely on a trustworthy calculation of the 
tunneling rates at present and in the early Universe. There is no reason why the Universe 
cannot be resting in a false vacuum which has a very long (on the cosmological scale) lifetime. 
We therefore disagree with the restrictions imposed by a number of authors |5], [7], |8[ on the 
allowed MSSM parameter space, which did not consider the corresponding tunneling rates. 
The calculation of the transition probability is more or less straightforward in the case of a 
single scalar but becomes extremely difficult for a potential that depends on several fields. 
Below, we address these difficulties and employ a new technique to determine the lifetime 
of the false vacuum in the case of the MSSM. 

We will see, in fact, that the SML minimum is effectively stable with respect to the 
transitions to the corresponding CCB minima for a substantial part of the allowed parameter 
space in the MSSM. We will also argue that, due to the specific nature of the CCB minima, 
they would not have been populated during the early stages of the evolution of the Universe, 
except for some small regions of parameters. On the other hand, the stability of the color 
and charge conserving vacuum on the time scales of order of the age of the Universe imposes 
important constraints on theoretical models and can provide guidance for future experimental 
searches. 
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The paper is organized as follows. In Section 2, we discuss general features of the SML 
vacuum decay in the MSSM. In Section 3, we consider the electroweak phase transition in 
the early Universe, in particular the issue of the 577(3) x SU(2) x t/(l)-symmetric vacuum 
stability before the transition. In sections 4 and 5, we study the zero-temperature tunneling 
rates numerically. The method used to compute the transition probability is described in 
Appendix A, while Appendix B contains an approximate description of tunneling in the limit 
of a very deep true minimum. 

2 Essential aspects of the MSSM vacuum stability 

We begin by considering a simplified version of the MSSM potential. As was emphasized in || 
^| , the third generation requires the most attention in connection with the issue of color and 
charge breaking, because the the CCB minima associated with the the large Yukawa coupling 
are the most dangerous^. We will see shortly that the tunneling rate into a CCB minimum 
is roughly proportional to exp(— c/y 2 ), where y is the corresponding Yukawa coupling and c 
is a constant. 

We begin by considering a model defined by the superpotential: 

W = yt L t R H 2 + f iH 1 H 2 , (1) 

where t L and t R denote the top quark superfields, and the Hi and H 2 are the MSSM Higgs 
bosons. At this point we ignore the leptons, lighter quarks, and the electrically charged 
Higgs components. The resulting scalar potential, including the soft SUSY breaking terms, 
is, at tree level, 

2 This statement was disputed in ||, where it was argued that in the limit of small Yukawa coupling 
both the cubic and the quartic terms are small and that only their relative values affect the depth of the 
CCB minimum. This is true, as long as one ignores the quartic D-terms which, in fact, become more 
important in the small Yukawa coupling limit. However, as was shown in Ref. fllpf , the height of the barrier 
separating the SML minimum from the CCB minimum is roughly proportional to l/y 2 ^, where y m in is the 
smallest Yukawa coupling associated with the fields that acquire non-zero vev in the CCB minimum. The 
corresponding tunneling rates are greatly suppressed for small y. This point is illustrated in the toy model 
discussed below (see also Ref. ||). Thus the CCB minima in which only the squarks of the third generation 
have non-zero vev play the most important role in our considerations. 
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where 



V = V 2 + V 3 + V 4 



(2) 



V 2 = m\Hf +m\Hf +2mlH° l Hl + mlP L +m 2 iR t 2 R (3) 

V 3 = -2AH° 2 tJ R -2^t L t R (4) 

V, = P L P H +? L H?+P a H? + V D (5) 

For the SU(3) x SU(2) x U(l) gauge group, the D-terms are 



V = ± 
D 8y 2 



gl{Hf - Hf - \? + \? f + gt(Hf - H? + P L f + ^ 3 2 (? " ? ) 2 



(6) 



Here the color indices are suppressed. We have absorbed the Yukawa coupling in equations 
(@-||) by the redefinition of the fields — > <j>/y and of the scalar potential V — > y 2 V. Also, 
all the fields are made real by a rephasing, and the complex phases are absorbed into the 
definitions of A and \x parameters. (There are strong experimental limits on such phases, 



which force A and /i to be nearly real; see Ref. |TTJ for reviews of these constraints.) 

In the i L = t R = hyperplane, equation (§) describes the usual MSSM Higgs potential. 
The constraint 



2 2 . / 2\2 , ( m l + m 2 
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mlml < (mi) 2 < ) (7) 

ensures the existence (first inequality) and stability (second inequality) of the minimum 
with the correct pattern of electroweak symmetry breaking. The latter inequality results 
from requiring that the quadratic term is positive definite along the flat directions of the 
quartic term V4. 

Much of our discussion will concentrate on the effects of the trilinear terms V3. If A 
and/or /x are large enough, the potential acquires an additional local, or global, minimum at 
some point outside the t L = i R = hyperplane. In this case, the electromagnetic U(l), color 



SU(3), as well as some other symmetries (e. g., the global £/(l)& aryon ) will be spontaneously 
broken by the non-zero vev of t L and i R . For example, the potential @) has, for appropriate 
values of A and fi, four degenerate CCB minima mapped onto each other by the following 
reflections: 
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The gauge SU(3) x SU(2) x Z7(l) symmetry is broken by the non-zero t L , £ fl and #2 
vacuum expectation values down to an 577(2) subgroup of the color SU(3). 

Evidently, for some otherwise reasonable values of the parameters, the potential may 
have a global CCB minimum. If this is the case, we would like to estimate the tunneling 
probability from the SML to the CCB minimum. 



The semi-classical calculation of the false vacuum decay width was done in Refs. [12, 13 



14 1 for the case of a single scalar field, 4>(x). For a recent review of tunneling we refer the 



reader to Ref. |[Lq| . The corresponding path integral is dominated by the field configuration 



4>{x) called the "bounce" and can be evaluated using the saddle point method |T4j . The 
bounce, being the stationary point of the Euclidean action, is the non-trivial solution of the 
corresponding Euler-Lagrange equation obeying certain boundary conditions. 
The transition probability per unit volume in the semi-classical limit WA 



is 



T/V = Ae- s[ ^ ]/h , (9) 

where S[4>] is the Euclidean action of the bounce, a classical solution to the variational 
equation 5S = 0. 

For the Universe to have decayed to the global minimum, the transition has to take 
place within a four-volume of size, roughly, t$, where to ~ 10 10 years is the age of the 
Universe. Taking the pre-exponential factor in ([5J) to be of order (100 GeV) 4 , one obtains 
(r/V)tp ~ 1 for S E [<f)]/h ~ 400. Therefore, a false vacuum whose decay rate is characterized 
by S E [4>]/h > 400 can safely be considered stable. 
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The presence of many scalar fields in the potential introduces a number of complications 
which will be addressed belowf| However, our immediate goal is to obtain a crude estimate 
of the false vacuum decay rate. We therefore make a further simplification and reduce the 
tunneling problem to that of a one-component case. 

Suppose the energy difference between the two minima AV = e 4 is small in comparison to 
the height of the barrier. Then the thin-wall approximation is appropriate and the Euclidean 
action of the bounce of size R is given by 

y 2 S E = -2n 2 ^-AV + 2n 2 R 3 f°° B sj2V(H) d0 (10) 

4 & SML 

where the appearance of the Yukawa coupling on the left-hand side of the equation is due 
to the fact that the fields have been scaled by a factor y in the beginning. 

The bounce (j)(x) corresponds to the extremum of the Euclidean action with respect to R, 
which is reached for the critical size of the "bubble" R c = 3Si/e 4 , where S\ = J ^2V(H) dcf). 
The corresponding action is 

97-7T 2 c ? 4 

Some comments are in order. First, we observe that, as was asserted above, the tunneling 
rate is very sensitive to the value of the Yukawa coupling. Therefore, the minima associated 
with the third generation of squarks are the most interesting. Second, it is well-known that 
the thin-wall approximation works well only for very small values of e. On the other hand, 
S B [<j>] < 400 in flnD corresponds to e/Sj /3 > [(27vr 2 /2y 2 )/400] 1 / 12 « 0.9, which is not the 
thin wall regime. This means that whenever the thin-wall limit is a good approximation, the 
transition we are interested in will not take place on the relevant time scale. The "dangerous" 
CCB minima lie, as a rule, outside the domain of validity of the thin-wall approximation. 

In Appendix B, we find an approximate representation of the bounce in the opposite limit, 
which we call a "thick- wall approximation" . Unfortunately, the phenomenologically accept- 
able values of the trilinear term can be approximated by neither thin-wall, nor thick-wall 
3 See also Refs. |,|l||l7|. 
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limiting expressions. For this reason, one must resort to a numerical analysis to determine 
the fate of the false SML vacuum in the presence of the lower lying CCB vacua. 

We will focus on the CCB minima associated with the scale of order the electroweak and 
the SUSY breaking scales, in contrast to other studies |L8], [RJ that discussed the possibility 
of a potential CCB minimum characterized by a planckian, or a GUT-scale scale vev. 

It was argued in Ref. |TB| that, if the MSSM is to be incorporated in some Grand Unified 
Theory (GUT) characterized by the scale M G , some new constraints should be imposed on 
the MSSM parameters to eliminate the possibility of a global CCB minimum developing at 
some scale Q, 1 TeV <C Q < M G . We note that symmetry restoration at a large energy scale 
is not required for the self-consistency of a spontaneously broken gauge theory. The latter is 
characterized by a symmetric action and an asymmetric vacuum. Cosmological data may, at 
least in principle, provide a test of whether the symmetric ground state existed in the Early 
Universe. However, the finite-temperature field theory describing an expanding universe has 
a different effective potential from that of the T = case. Most of the flat directions of 
the tree-level potential are lifted by terms of order SUSY breaking scale []. Therefore, the 



CCB minima of the kind dealt with in Ref. [18| may not be present in the Early Universe. 



Furthermore, as was pointed out in Ref. ||19|| , the existence of a global CCB minimum at such 
a large value of the vev is irrelevant for the low-energy physics because the color and charge 
conserving vacuum is effectively stable with respect to tunneling into the CCB minimum of 
that kind. 

This is an example of a phenomenon which is similar to the usual (perturbative) de- 
coupling. Tunneling into a "very deep" vacuum occurs at a rate which is independent of 
the depth of such a minimum and is determined only by the magnitude of the field and 
the steepness of the potential at some well defined point, the "escape point", (fi e (c. /. the 
discussion in our Appendix B). The shape of the effective potential at the scale Q > ^ e 
has no effect on the transition probability, and therefore the low-energy physics is indepen- 
dent of the physics at Q 3> <fi e - This is analogous to decoupling in perturbation theory, even 



4 In addition to finite-temperature and tree-level breaking of supersymmetry, there are SUSY breaking 
terms associated with the metastability of the false vacuum [§3 . 
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though the perturbative decoupling theorems may not apply to non-perturbative effects such 
as tunneling. 

In practice, this decoupling allows us to treat the ultra-deep minima on the same footing 
as the CCB minima of the tree-level potential. Also, since all the relevant dimensionful 
quantities, including e , are of order 100 GeV to 1 TeV, the radiative corrections to the 
effective potential have a very small effect on the tunneling rates. Therefore, it is well- 
justified to use a tree-level potential for evaluating the stability of the SML vacuum and 
determining the allowed regions of parameters in the MSSM. In a future study, we plan 
to further refine our present results by taking into account the radiative corrections to the 
effective potential. 

Closely related is the issue of the directions in the scalar sector of the MSSM along which 
the tree-level potential appears to be unbounded from below (UFB). These directions are 
chosen to zero the quartic terms in the tree-level potential. Usually, the one-loop radiative 
corrections rectify the situation by introducing positive definite quartic terms of the type 
STrM A \og(M 2 /Q 2 ). If this is the case, the full effective potential turns out to have only a 
very deep CCB minimum and is not unbounded from below. The latter may be separated 
from any other vacuum by a high enough barrier to make the presence of such a CCB 
minimum irrelevant. To determine whether a certain region of the parameter space must be 
excluded, one must again examine the corresponding tunneling rates. However, since in any 
case the tunneling rate is determined by the shape of the potential at the vev's of order a few 
TeV, one can treat the UFB directions as if they were leading to a very deep minimum. We 
stress that this is true, regardless of whether the given direction is a UFB direction of the 
exact effective potential, or only of its finite-order (tree-level, one-loop, etc.) approximation. 

3 CCB minima in the early Universe 

If the effective potential has more than one minimum, then the determination of the physical 
vacuum requires consideration of the evolution of the Universe. For the MSSM, the analysis 
is complicated by the presence of flat directions (see, e. g., Ref. f2l|), lifted only by terms of 
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order the SUSY breaking parameters, along which the scalar field may acquire a large vev 
and fall out of thermal equilibrium in the early Universe. We leave this question to future 
study and assume that the electroweak phase transition proceeds in the usual manner, from 
an SU(3) x SU(2) x [/(l)-symmetric to a broken phase. 

This assumption is well justified for the types of CCB minima we consider. It is well- 
known that in inflationary models the large fluctuations of the scalar fields may populate 
some color and charge breaking minimuma. This is true of both local and global minima. If 
the vacuum expectation value of the scalar field in that minimum is large in comparison to 
the reheating temperature, the Universe may "freeze" in that minimum. However, the CCB 
minima we consider, unlike those of Refs. [IR ITTJ], have vev's of the order of the electroweak 
scale and usually disappear at temperatures of order 1 TeV, except for the models with 
non-trivial (and atypical) symmetry restoration pattern like that of Ref. 0. Since most 
inflationary models predict much higher reheating temperatures, one can assume that at 
T ~ 1 TeV the SU(3) x SU{2) x U(l) is unbroken. 

The main effect of the temperature-dependent corrections to the effective potential is 
contributions of order T 2 to the quadratic terms. The trilinear terms also receive some 
corrections, linear in T. 

The depth of the CCB minima depend on the relative values of the squark mass terms 
(mo) on one hand, and A or /i on the other. At finite temperature, positive mass-squared 
terms proportional to T 2 appear in the effective potential and lead, at some critical temper- 
ature, T c ~ 100 GeV, to the disappearance of the SML minimum. Since A, \x and tuq are 
allowed to be large in comparison to T c , it is possible that as the Universe cooled the nega- 
tive energy CCB minima formed at some T CCB > T c , before the Universe was cold enough 
to undergo the electroweak phase transition to the SML minimum. This would allow for 
the possibility of a transition from an SU(3) x SU(2) x £/(l)-symmetric to the color and 
charge breaking phase. One then must consider whether the transition to the CCB minimum 
actually occurred, and, if so, what happened subsequently. 

Then, a priori there are three possibilities: (i) the Universe may stay in the symmetric 
minimum until the temperature reaches T c and the usual electroweak transition takes place; 
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(ii) the Universe may go to the CCB minimum and freeze there; and (hi) the transition from 
the symmetric to the SML minimum may occur in two stages: first the transition from the 
unbroken to the CCB phase takes place and only then the Universe may tunnel to the SML 
vacuum. Clearly, the second possibility is excluded empirically. 

Possibility (iii) is intriguing. The idea of a multi-stage phase transition, in the course 
of which the gauge group might change a number of times finally arriving at the Standard 
Model group, is not new. An evolution of this kind could have important implications for 



magnetic monopoles [p2| , charge asymmetry of the Universe and magnetic field generation 
2~3fl , baryogenesis p4j, etc. The sufficient (though not necessary) conditions for the first 
part of this scenario, the transition from the symmetric to the CCB phase, will be derived 
below. However, we have not found a model in which a second-order phase transition from 
the CCB minimum to the SML minimum would subsequently take place. It appears to be a 
generic feature that even when the CCB minima can be populated via a second-order phase 
transition, they are separated from the SML minimum by a thick barrier. If the Universe is 
stuck in the false CCB vacuum (which, in this scenario, must have a higher vacuum energy 
than the SML phase), this would trigger a new stage of inflation which could only be ended 
by a first order phase transition characterized by a relatively low tunneling rate. It was 
shown in Ref. |25| that this sort of transition proceeds via bubble nucleation at a rate 
which is too slow to catch up with the expanding volume. The bubbles of the true vacuum 
would neither collide, nor percolate, preventing the Universe from reheating. Therefore, the 
two-stage electroweak phase transition could not have taken place in the early Universe. 

Let us consider the possibility that a CCB minimum exists at some temperature T > T c . 
At that temperature, the mass matrix of the third generation squarks (ignoring the rest of 
the squarks and sleptons) in the 577(2) x [/(l)-symmetric phase (Hi = H 2 = 0) is of the 
form: 



ml +c r T 2 

,2 



(12) 

Suppose that it has at least one negative eigenvalue, which makes the SU(2) x £7(1)- 



mf +c R T 2 

t R 



10 



symmetric vacuum unstable. Since T > T c , the SML minimum does not yet exist. The 
decay of the unstable vacuum will result in the creation of the CCB condensate by a second- 
order phase transition. At the same time, we assume that all the necessary conditions have 
been applied to constrain the masses- squared of squarks to be positive (and large enough) 
in the SML minimum at zero temperature, where the mass matrix is 



2 _i_ oto2 

.-i l l ) 7 II I I 1 

(13) 



mf + mf AH 2 + fiH 
AH 2 + fiH 1 mf + m, 

R 




where the mt is the top quark mass. 

We would like to see whether there is a region of parameter space in which the matrix 



(|13|) has only positive eigenvalues, while the matrix (12) has a negative eigenvalue. It is 
easy to see that this is only possible if either c L T c 2 < mf, or c H T c 2 < mf. The shaded region 
shown in Figure [I] corresponds to the additional domain of parameters which can be ruled 
out by requiring stability of the SU(2) x U(l) symmetric minimum above the electroweak 
transition temperature. This domain comprises two regions in which mf and mf differ in 

t L R 

sign (Figure [l]). The hyperbola (mf +mf)(mf +mf) = (AH 2 + fiHi) 2 , where the values of 

L R 

Hi and H 2 are computed at the SML minimum, outlines the domain of positive determinant 
of the matrix in equation (p!3|). If one requires that the squark masses be greater than 45 
GeV, in accordance with the current experimental limits, it would further reduce the area 
of the corresponding (shaded) regions in Figure |l|. 

The coefficients c L and cr are of order 1; their exact values depend on the spectrum of 



the MSSM and have been computed, e. in Ref. [26|. For new constraints to arise from 
the requirement of symmetric vacuum stability at T > T c , the inequality T c ~ 100 GeV < 
m t/\/ c L,R mus t be satisfied. 

We note that the high-temperature expansion |2^J is not expected to be accurate for 
T ~ T c . Therefore, the accuracy to which one can determine the boundaries of the shaded 
regions in Fig. [j] is limited by one's inability to determine the effective potential accurately 
for T ~ T c due to the limitations of the theoretical framework |27| . It is clear, however, that 



the size of the shaded regions in Fig. [L] is rather small, so their exclusion cannot cause an 
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Figure 1: Region of parameters which can be ruled out by requiring stability of the SU(2) x 
£7(1) symmetric minimum above the electroweak transition temperature. 

appreciable reduction in the MSSM parameter space. For the generic values of parameters 
outside the shaded regions, the second-order phase transition into the CCB minimum does 
not take place. 

We have demonstrated that the symmetric vacuum of the MSSM is generally stable with 
respect to second-order phase transitions to a CCB minimum at T > T c . However, the 
question of a first-order transition at finite temperature remains open. The probability of 



tunneling in the high-temperature limit is suppressed by the factor |28 

r/Voce~ s ^ hT 



(14) 



where 5*3 [(f)] is the three-dimensional action of the d = 3 bounce. The time allowed for the 
transition is roughly t T = m pl /T 2 , the age of the Universe when the temperature equals T, 
which means that S^[(j)\/hT must be less than about 45 for the transition to take place. 

Thus, it appears most likely that if T c is not too small (of order 100 GeV, as is generally 
believed), then the second-order transition to the CCB minimum of the type discussed 
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above is not likely, and the Universe is driven towards the color and charge conserving SML 
minimum of the scalar potential. While, admittedly, this is not a rigorous theorem, the 
color and charge conserving minimum appears to be favored by the thermal evolution of the 
Universe. We leave the detailed investigation of this issue for future work. 

Our next question is what happens after the Universe cools down in the SML minimum 
of the potential: will the false vacuum be effectively stable at T « even in the presence of 
some deeper CCB minima, or not. 

4 Numerical analysis of tunneling rates at zero tem- 
perature 

To make a definitive determination of whether the SML vacuum is stable with respect to 
decay into a lower CCB minimum, one has to compute the tunneling probability numerically. 
Analytic computation is usually feasible only in the thin wall limit. 

In the semi-classical limit, the zero-temperature tunneling probability p| |T^, [TJ|, [H], 17 1 
per unit volume is 



r/v = c G 



(x)\ 2 d 4 x 



1 N / 2 



2nh 



det'[-dl + U»($)] 



det[-02 + U" 



-1/2 



x(l + 0(ft)) (15) 



where S[4>] is the Euclidean action of the bounce, a solution to the variational equation 
SS = 0, det' stands for the determinant with all the zero eigenvalues omitted, N is the 
number of Goldstone zero modes and the C G is the group-theoretical coefficient JTBL IT7J. 



Suppose the scalar potential U(<pi, <f> n ) has a local minimum at 



l,2,...,n, 
l,2,...,n; 



as well as at least one additional (local or global) minimum at <pi = <j>% i 
[/((/>*) < U ((/>*). Then the bounce 4>(x) = (<f)i(x), 4> n (x)) is a non-trivial 0(4)-symmetric 
solution 0(r), r = y/x^, of the system of Euler-Lagrange equations: 



d - - 
-^[7(01, . ..,</>„) 



(16) 



with the following boundary conditions: 
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' (d/dr)&(r)\ r=0 = 

< _ (17) 
4 0i(oo) = (f) f 

In the case of a potential that depends on a single scalar field, one can solve equation 
(PI) with the boundary condition ( |TTD numerically. The straightforward technique is to 
assume some value for the unknown quantity (f)(0) = <ft e , the so called "escape point". Then 
one can integrate equation (|T6| ) numerically and vary the escape point until the proper 
limit 0j(oo) = <fif is reached. Here one uses the fact that if for some value <p e = <p e _ the 
corresponding 4>i(oo) > , and for some other value e = 0j(oo) < <f)f , then the true 
escape point must lie somewhere between 01 and <p e + . This strategy fails when one has to deal 
with a potential that depends on several scalar fields. The peculiarity of one- dimensional 
topology no longer allows one to find the true escape point as the compromise between the 
"undershoot" and the "overshoot" Q. It is a generic property of the bounce that it is a saddle 
point of the Euclidean action [|lj| , and not a minimum, and thus small changes in the initial 
conditions result in large changes in the form of the solution. Therefore, it is impossible to 
find the bounce numerically in the case of a mult i- component field using the procedure just 
mentioned. 

This difficulty was realized by the authors of ||, who proposed an iterative procedure 
to look for the bounce as a special point of the discretized action on a lattice. Since the 
solution in question is not a minimum, but a saddle point, they tried to minimize the action 
with respect to random variations, while maximizing the same action with respect to scaling 
r — > Ar. We find that the iterative procedure of this kind is ill-defined and cannot have a 
meaningful limit. First, it is impossible to separate the variations corresponding to scaling 
from those orthogonal to scaling in a practical numerical simulation. Second, although it 
is true that the bounce maximizes the action with respect to scaling, it is easy to see that 
the generator of such a variation cannot be the eigenvector of the second variation operator 
corresponding to the negative eigenvalue. 

5 Even in the case of a single scalar field, the shooting method may turn out to be ineffective for a potential 
with sufficiently degenerate minima. This is because one may be required to specify the trial value for the 
escape point with an exceedingly high precision to compute the action of the bounce to a given accuracy. 
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Instead, we use a new method proposed in || to find the solution of fll6D ,fll7D. The 
idea is to turn the saddle point of a Euclidean action into a true minimum by adding to 
it some auxiliary terms which vanish for <ft = <j)(x). The resulting "improved action" will 
have a minimum at the point corresponding to the desired solution, which can now be 
found by minimizing the discretized version of improved action on the lattice. Details of the 
application of this method are given in Appendix A. 

5 The MSSM vacuum stability at zero temperature 

For the reasons explained above, we consider the MSSM potential with the third generation 
only. We also neglect all of the trilinear couplings except those proportional to the largest 
Yukawa couplings, y t and y^. This justifies dropping the squarks of the first two generations. 
Also, although we allow the values of tan/? from ~ 1 to ~ 60, the most stringent constraints 
come from the small tan/3 region, where the top Yukawa coupling is larger. As was pointed 
out in Ref. ||, the noteworthy CCB minima can only develop in the directions along 
which one and only one of the trilinear terms is non-zero. Thus, it is the size of the largest 
Yukawa coupling that affects the lifetime of the false vacuum the most. We therefore consider 



the following subset of the MSSM scalars: $ = H$, iff, H% , Q L , t R , b R }, where Q L = 



corresponds to the following scalar potential (with the soft SUSY breaking terms included) 
at tree level: 



The superpotential 



W = y t Q L H 2 t R + y b Q L H l b R - y,E x E 2 



(18) 



v = v 2 , H + v 2 + v 3 + V A , 



(19) 



where V 2 jj comprises the terms quadratic in the Higgs fields, 




(20) 
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V 3 = - yt A t {Hlt L -HfbJ R - yb A b {H^t L -Hlb L )b R 
-y t n(H°ytj R - y bf i(H+)%b R 

-ytii(HD*b L i B - y bf i(H°y~b L ~b R + h. c, (21) 

and V4 is comprised of the D-terms, as well as the terms of the form yl b <^\^ where <pi G 

We now come to the issue of radiative corrections. It has been argued (see, e. g., 0) that 
if the CCB minimum occurs at a scale which is within one or two orders of magnitude from 
the electroweak scale, then the (logarithmic) radiative corrections are too small to cause a 
noticeable distortion in the shape of the effective potential, and therefore can be ignored. 
While this conclusion is correct, we would like to emphasize that, as was discussed in Section 
2, it is not so much the position of the CCB minimum, as the size of the bounce and the 
value of the escape point that determine the scale to which the bounce will "probe" the 
effective potential. The shape of the effective potential beyond that scale has no effect on 
the tunneling rate. In our numerical calculations, a typical size of the bounce for realistic 
values of of the MSSM parameters is of order (1 — 0.1 TeV)" 1 (except for deliberately fine- 
tuned cases of highly degenerate minima in which the bounce blows up to the size of its 
thin- wall limit, R ~ 1/e, and the tunneling becomes highly improbable). This justifies post 
factum the neglect of the one-loop contribution in the effective potential. As long as one 
never encounters very small (in length units) bounces, one can be certain that only the shape 
of the potential around the electroweak scale is relevant for tunneling. 

We search the MSSM parameter space by generating randomly the values of the param- 
eters that enter in the scalar potential (|T9|) . Then the minima are found numerically and 
the tunneling probability is computed using the method of Appendix A. 

The results of the numerical analyses are plotted in Figure 0, where the action of the 
bounce S is shown as a function of the relative depth of the CCB minimum with respect to 
the SML minimum, while the other MSSM parameters are varied randomly. For each point 
plotted, the global minimum of the potential is not color and charge conserving. Nevertheless, 
one observes that S takes values on both sides of the critical value, S = 400. Therefore, for 
some values of parameters, namely those for which S > 400, the false SML vacuum is stable 
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Figure 2: The action of the bounce, S, does not depend sensitively on the depth of the CCB 
minimum, AV^, except in the "thin- wall" limit (small Ay). 

on a time scale large compared to the age of the Universe. Such color and charge breaking 
minima are "safe", and the corresponding parameters are allowed. In contrast, those points 
that correspond to S < 400 are ruled out by the mere existence of the world as we know 
it. Another lesson one learns from Figure ^] is that the tunneling rate is not very sensitive 
to the depth of the CCB minimum, except in the limit of nearly degenerate minima, the so 
called "thin- wall" limit. This is to be expected, as was discussed in Section 2. 

The domain of stability of the false SML vacuum with respect to tunneling is delineated 
by stars in Figure [| The lighter top squark in the presence of the large trilinear couplings 
forces the barrier, which keeps the system in the metastable vacuum, to be thinner and 
lower. That results in a higher likelihood of tunneling. The points labeled by boxes fall into 
the domain that is excluded by the existence of our (color and charge conserving) Universe. 
On the other hand, if both left-handed and right-handed stops are heavy, and if the trilinear 
terms are small, then the false vacuum is stable and the presence of a global CCB minimum 
is irrelevant. We note that a number of models favor the left lower corner of the plot in 
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Figure |3], where only a few "dangerous" CCB minima occur. 

It would be useful to derive some empirical algebraic constraints to distinguish between 
the allowed and excluded domains of parameter space based on the numerical results. For 
instance, it has been argued || that the inequality 

4 2 + 3 / /<3(m? L +m?J (22) 

is a reasonably good condition for excluding a global CCB minimum. This inequality would 
force one to be above the dotted line in Figure [|. This is in good agreement with our 
numerical results. However, if one takes into account the tunneling rates, the constraint ( p2|) 
is relaxed significantly. The empirical inequality which should replace (|22"D, 

A 2 t +3fi 2 < 7.5 (m\ +mf K ), (23) 

is depicted by a thick dashed line in Figure §. The simple condition (^3|), should be applied 
with caution because, strictly speaking, it is neither necessary nor sufficient. It is an ap- 
proximate empirical inequality that may be useful for a crude determination of whether the 
CCB minima are "dangerous". 

For the phenomeno logically attractive values of < 2 TeV, \A\ < 4 TeV, it is generally 
true that "the larger the trilinear coupling, the more dangerous is the corresponding CCB 
minimum" . However, it is instructive to examine what happens to the tunneling probability 
in the limit of very large fi and A t (and large enough squark mass terms to ensure the 
existence of the SML minimum). In that limit, as the CCB minimum moves away from the 
SML minimum, the barrier separating the two becomes thicker, and the false vacuum should 
become more stable. This is, in fact, what happens. The set of point in Figure [5] includes 
those points (located in the lower left corner) shown in Figures [|, |3] and [|. In addition, Figure 
H displays the points corresponding to some very large values of A t and /i. As expected, the 
tunneling probability diminishes for very large values of A t and /i, and and m$ . 

To summarize, if the global CCB minimum is nearly degenerate with the local SML 
minimum (thin- wall limit), then the tunneling probability is extremely small. As the trilinear 
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Figure 3: The domains of stability (stars) and instability (boxes) of the false SML vacuum 
with respect to tunneling into the global CCB minimum. Light top squark and large trilinear 
couplings generally correspond to a lower and thinner barrier and, thus, higher probability 
of tunneling. 
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Figure 4: Each point represents the set of the MSSM parameters for which the global 
minimum of the potential is color and charge breaking. The stars correspond to the SML 
false vacua whose lifetime is large compared to the age of the Universe. The boxes indicate 
those points in the parameter space for which the false SML vacuum should have decayed 
via quantum tunneling. The dotted line represents the empirical criterion for the absence of 
the global CCB minima: A 2 + 3u 2 < 3M 2 , where M 2 = m 2 + m 2 . Taking into account the 

tunneling rates relaxes this constraint to, roughly, A\ + 3/i 2 < 7.5M 2 , shown as the dashed 
line. The scale is logarithmic. 
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Figure 5: Tunneling probability for unphysically large values of A t and fi. As the CCB 
minimum moves farther away, it becomes "less dangerous". As before, the stars mark the 
points with S > 400, while the boxes depict those with S < 400. 



couplings increase, the false vacuum decay rate increases because the escape point of the 
bounce moves out of the flat vicinity of the global minimum into the region in which the 
gradient of the potential is significant. However, a further increase in the size of the trilinear 
couplings, as well as the consequent increase in the squark mass terms, makes the barrier 
thicker and pushes the escape point away from the SML minimum. This eventually causes 
a decrease in the tunneling rate. In accordance with one's intuition, the low-energy physics 
is unaffected by the physics at the very high energy scales. 
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6 Conclusion 



The color and charge conserving minimum may not be the global minimum of the MSSM 
potential. It is possible that the Universe rests in a false vacuum whose lifetime is large in 
comparison to the present age of the Universe. Under fairly general conditions, the SML 
vacuum may be favored by the thermal evolution of the Universe, even if it does not represent 
the global minimum. 

The existence of the CCB minima of the scalar potential results in some important 
constraints on models with low-energy sypersymmetry. However, the commonly imposed 
(see, e.g., Refs. |3(| |3l] ) requirement that the SML minimum be global is too strong and 
may overconstrain the theory. In fact, for a large portion of the parameter space the presence 
of the global CCB minimum is irrelevant because the time required for the Universe to relax 
to its lowest energy state may exceed its present age. The basic reason for this is that the 
quantum tunneling is a non-perturbative phenomenon that is naturally associated with the 
energy scale that is exponentially smaller (suppressed by a factor exp{— S}) than the typical 
scale in the theory. 

We have computed numerically the SML false vacuum decay rates for a variety of values 
of the MSSM parameters. Our results indicate clearly that the MSSM vacuum stability 
with respect to tunneling into a CCB minimum or a "UFB direction", imposes important 
constraints on models with low-energy supersymmetry. 
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Appendix A 



The numerical algorithm for computing the tunneling probability comprises several steps. 
First, the positions of the minima, 0* and <\>* of the potential U((f>) are determined numerically 
for a given set of values of the MSSM parameters. If there is no minimum below the SML 
vacuum, a new set of parameters is chosen. We define the Euclidean action on the L-point 
lattice: 



S[d>] = T[<j>) + V[<f>} (24) 
T[0] = 2vr 2 A 4 EK + i-^ E^S (25) 

m=l \i=l *\ a m+l - a m) A J 

V[<f>) = 2n 2 A'J2(d m+1 -d m )d 3 m U(^ l \---A ( r y ) (26) 

m=l 

where d m is the position of the m's site of the lattice in dimensionless units, and A is 
the length parameter that determines the overall scale and is chosen so as to optimize the 
computation. 

The next step is to define the improved action || , for which the bounce is the minimum. 
We do that by adding to the action some auxiliary terms which (i) vanish as <fi{x) approaches 
the bounce, and (ii) make the bounce 0(r) a minimum of the improved action ||: 

~S[<p] = S\<P} + \ |T[0] + 2^[0]| 1/2 (27) 

where A is some arbitrary (dimensionless) Lagrange multiplier. 

The effect of adding the auxiliary terms is shown qualitatively in Fig. |], where two 
projections of the saddle point of S[4>] are depicted symbolically. 

The bounce 4>(x) is always a minimum of the improved action S[<p]. However, unless A 
is chosen to be large enough, the difference S[(j>] — S[(p] may appear to be too small in the 
vicinity of the bounce. In this case, it is possible that the small perturbations will take one 
over the barrier (Fig. ^ towards the trivial solution of zero action 4>(x) = $ . Thus, it is 
crucial to take A large enough for the numerical simulation to succeed. 
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Figure 6: The improved action S[(f>] (dashed curve) has a minimum that coincides with the 
saddle point of the actual Euclidean action S[(p] (solid curve), as shown in two projections. 
The origin on the right corresponds to the trivial solution <f>(x) = $ for which S[(j>] = S[(f)} = 
0. The size of the "gap" between the two curves is a function of A. 
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Initially, one may choose 



Am) _ j $, m<L/2 

~ { <t>l m > L/2, (28) 

where 0° is a zero of V (<j)). Alternatively, one may start from a different profile for the 

bounce, e. g. the thick-wall ansatz described below in the Appendix B. 

Then we find the optimized value for the overall scale A by maximizing the action with 
respect to A. After that, we allow random variations of the values of the scalar field at 
each lattice site to minimize the improved action (p7[). The iterations stop when further 
variations do not lead to a reduction in the improved action. As an independent criterion 
of the quality of the fit, we require that T[0]/(— 2 V [</>]) = 1 ± 0.01 (see, e. g., Ref. and 
references therein for a discussion of this identity). 

As was explained in Sections 2 and 5, tunneling into a "UFB direction" is equivalent 
to a transition into a very deep CCB minimum. The details of the effective potential at 
scales much larger than e do not affect the decay probability. In particular, it is irrelevant 
whether the alleged "UFB direction" is heading towards a very deep minimum, or minus 
infinity. This allows one to introduce an effective cutoff to stop the run-away fields from 
going to infinity. In practice, we did not allowed the value of <ft to run beyond 10 TeV. 
Therefore, the "UFB directions" were treated as if they lead to a CCB minimum with a 10 
TeV vev. In each particular case this procedure is justified post factum by ensuring that the 
value of the escape point is small compared to a 10 TeV cutoff. 

Appendix B 

In the limit of nearly degenerate minima separated by a high barrier, the so-called thin- 



wall limit, the bounce can be approximated by a smoothed out step function fll3[ . This 
approximation proved very useful in estimating the tunneling rates and was used in our 
analysis of the toy model in Section 2. In practice, however, one rarely encounters a situation 
in which the thin- wall approximation is in good agreement with the numerical results (c.f. 
Ref. ||). The tunneling rates are usually very small in the thin- wall limit, and therefore 
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in many physically interesting models the first-order phase transition takes place when the 
energy difference between the true and the false vacuum is not small in comparison to the 
height of the barrier. Then it is necessary to go beyond the thin-wall approximation, which 
is usually done by means of a numerical calculation. 

If, however, the energy difference AV between the two vacua is larger than the height of 
the barrier, one can find a simple approximation to the bounce in what we call a "thick-wall" 
limit. As AV increases, the so-called "escape point" , the value of the field in the center of the 
bounce <p e = 0(0), moves away from the minimum of the potential. Then equation ( ]T6|) can 
be linearized in the vicinity of the escape point because dU(<p)/d(f) is a constant independent 
of 4> for <p ~ 0(0). Therefore, in the vicinity of the center of the 0(4)-symmetric bounce 



0"(r) + f 4>'{r) = -a 

(29) 

0'(O) = 0, 



where the constant a 



dU i 



. This is a linear equation whose only solution satisfying the 



boundary condition is of the form 

0( r ) = _ V + e (30) 



Outside the small neighborhood of the origin, the approximation fl29|) is not valid and 
(f){r) falls off exponentially, just as in the thin- wall limit ||13|| . The equation for the bounce 
0(r) for large r (and, therefore, small </>) becomes 

' 0"( r ) = dU{(j))/d<j) « m 2 2 

(31) 

, 0(oo) = 0, 

where m is the mass term for the potential U(4>) = (m 2 /2)0 2 + .... The solution of the 
equation (|31| ) is 

0( r ) = Ce~ mr , (32) 



2(3 



where C is an arbitrary constant. We now sew the approximate solution for the bounce 



from the two asymptotics, fl30|) and (p2[) , at some point r = R. The values of C and R 
are determined by requiring continuity and differentiability of the solution at r = R. The 
resulting ansatz is 

e - (a/8)r 2 , r < R 
0(r) = < (33) 



aR^—m(r—R) r ^> R 
4m ' — ' 



where 



R = —[yfi + 8m 2 e /a- 1], (34) 
and the value of e is an unknown parameter which can be found either from requiring that 

d(f) e d(fi e ' ^ 

or, equivalently, from solving the equation 

T = -2V, (36) 



where 



T = 2tt 2 /°° r 3 dr - ( ^^-) and 7 = 2vr 2 [°° r 3 dr UU(r)) (37) 
J o 2 \ dr J Jo 

are functions of e . 

The representation of the bounce in the thick-wall limit described above is approximate 
and is not very useful in application to the MSSM. However, it exhibits the essential features 
of tunneling in this limiting case. We also found it convenient to use this approximate solution 
as an initial profile for the bounce in the numerical procedure described in Appendix A. 
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